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ABSTRACT 

Within the framework of Newtonian magneto-sohd-mechanics, relying on 
equations appropriate for a perfectly conducting elastic continuous medium 
threaded by a uniform magnetic field, the asteroseismic model of a neutron star 
undergoing axisymmetric global torsional nodeless vibrations under the combined 
action of Hooke's elastic and Lorentz magnetic forces is considered with empha- 
sis on a toroidal Alfven mode of differentially rotational vibrations about the 
dipole magnetic moment axis of the star. The obtained spectral equation for 
frequency is applied to £-pole identification of quasi-periodic oscillations (QPOs) 
of X-ray flux during the giant flares of SGR 1806-20 and SGR 1900+14. Our 
calculations suggest that detected QPOs can be consistently interpreted, within 
the framework of this model, as produced by global torsional nodeless vibrations 
of quaking magnetar if they are considered to be restored by the joint action of 
bulk forces of shear elastic and magnetic field stresses. 

1. Introduction 

The recent detection of quasi-periodic oscillations (QPOs) on the light-curve tails of 
X-ray flaring SGR 1806-20 and SGR 1900+14 (Israel et al 2005; Watts & Strohmayer 2006) 
that have been ascribed to torsional seismic vibrations of quaking magnetar raises sev- 
eral questions of quite general interest for the asteroseismology of degenerate solid stars. 
Namely, whether these differentially rotational oscillations are predominately of an elastic 
nature, that is, restored by Hooke's force of shear solid-mechanical stresses or they should 
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be thought of as the toroidal magnetic Alfven mode of axisymmetric vibrations about dipole 
magnetic moment axis of the star and restored by the Lorentz force of magnetic field stresses 
(Glampedakis, Samuelsson & Andersson 2006). Also, it remains questionable whether these 
vibrations are of global character, that is, excited in the entire volume or they can be ex- 
plained as locked in the peripheral finite-depth crustal region of a neutron star. These and 
related issues are currently the subject of intense theoretical investigations by classical the- 
ory methods of material continua (Piro 2005; Levin 2007, Lee 2008; Bastrukov et al, 2007a; 
2008a) and by general relativity methods (Sotani, Kokkotas & Stergioulas 2008; Samuelssen 
& Andersson 2007). 

In this article, continuing the above investigation, we study an asteroseismic model of 
a neutron star with uniform internal and dipolar external magnetic field undergoing quake- 
induced global differentially rotational, torsional, vibrations about the dipole magnetic mo- 
ment axis of the star under the joint action of the elastic Hooke's and the magnetic Lorentz 
forces. In so doing we confine our attention to the regime of extremely long wavelength 
differentially rotational fluctuations of material displacements about motionless stationary 
state which are insensitive to core-crust compositional stratification of quaking neutron star. 
The characteristic feature of this regime is that conducting and highly robust to compres- 
sional distortions of solid-state material, of both core and crust, coupled by Maxwell stresses 
of uniform fossil magnetic field frozen in the star, sets in coherent axisymmetric differentially 
rotational vibrations with the nodeless toroidal field of material displacements identical to 
that for global torsional vibrations of a spherical mass of an elastic continuous medium ca- 
pable of transmitting perturbation by transverse shear elastic waves generic to solid state 
of condensed matter, not a liquid on^. In the context of the above QPOs problem, such 
vibrations have been analyzed in some details in our recent works (Bastrukov et al 2007a, 
2007b, 2008a), but the driving force was assumed to be of only elastic nature, that is, owing 
its origin to fluctuations in shear elastic stresses. The focus of this paper is placed, therefore, 
on the toroidal Alfven vibrational mode in which magnetic fleld and fleld of material dis- 
placement undergo coupled fluctuations restored by Lorentz force with Amper's form of the 



^It is commonly agreed that the average speed of transverse wave of elastic shear in the crustal matter is 
Ct = \/ m/ P ^ 10^ cm s""'^ which holds for both the outer crust with the density 10^ < P < 10^"'^ g cm"-^ and 
the inner crust, whose density is ranged between the neutron drip line density p ^ 10^^ g cm"'^ and density 
of overwhelming neutronization p ^ 10^"^ g cm~^ (e-g-i Blaes et al 1990). Accordingly, the shear modulus of 
the crustal matter is ranged from the surface to the crust-core interface as 10^"* < p< 10^^ dyn cm . The 
much denser core matter compressed by self-gravity to nearly normal nuclear density p = 2.8 • 10^* g cm~^ 
is much harder to elastic shear distortions and its shear modulus is estimated as /i ~ 10'^^ dyn cm~^ (e.g., 
Owen 2005, Bastrukov et al 2008b and references therein); in view of this the Fermi-degenerate neutron 
matter of atomic nuclei and neutron star cores is designated as an elastic Fermi-solid capable of transmitting 
shear elastic wave with the speed Cj = \fpjp ~ 0.2 c where c being the speed of light. 



-3- 



conduction current density. Before proceeding to the details of calculations it worth noting 
that the liquid star model with axisymmetric poloidal uniform magnetic field inside and 
dipolar outside has been the subject of serious investigations in the past in the context of 
magnetic variables (e.g., Schwarzschild 1949; Chandrasekhar & Fermi 1953; Ledoux & Wal- 
raven 1958). The situation may be quite different for ultra strong internal magnetic fields 
frozen-in the super dense matter of the end products of stellar evolution, white dwarfs and 
neutron stars. In neutron stars self-gravity is counterbalanced by the degeneracy pressure 
of relativistic electrons in the crust, whose highly conducting matter is capable of sustaining 
persistent current-carrying flows and by degeneracy pressure of non-relativistic neutrons in 
the cores of neutron stars whose material can be in the state of paramagnetic magnetization 
caused by Pauli's mechanism of field-induced alignment of spin magnetic moment of neutrons 
along the frozen-in the star fossil magnetic field (e.g., Bastrukov et al; 2002a; 2002b). Of 
course, there is no compelling evidence that this is the case and there is no certainly general 
agreement how the magnetic fields of pulsars and magnctars arc produced. Over the years a 
model of a neutron star with uniform internal and dipolar external field, pictured in Fig.l, 
has been invoked to discussions of evolution of magnetic dipole fields of isolated radio pulsars 
(Flowers & Ruderman 1977), whose surface dipolar magnetic fields are found to be highly 
stable to spontaneous decay (Bhattacharya & van den Heuvel 1991, Chanmugam 1994), as 
well as of origin and evolution of ultra strong magnetic fields of magnetars (Braithwaitc & 
Spruit 2006; Geppert & Rheinhardt 2006; Spruit 2008), therefore the study of axisymmetric 
torsional vibrations within the framework of a neutron star model with uniform internal mag- 
netic field seems amply justified. This interesting in its own right model describes situation 
when physically meaningful analytical solution of the eigenfrequency problem can be found 
and discloses mathematical difficulties one must confront when computing the frequency of 
Alfven vibrational modes in the neutron star models with non-homogeneous internal and 
dipolar external magnetic fields. An analytic example of such latter field is discussed in 
(Roberts 1981; Geppert & Rheinhardt 2006). 

The paper is organized as follows. In Sec. 2, the exact solution is given for the eigen- 
frequency problem of global torsional nodeless vibrations of a solid star with a uniform 
field inside and dipole outside driven solely by Lorentz restoring force. The obtained fre- 
quency spectrum of toroidal and poloidal magnetic Alfven vibrational modes in such a star is 
then quantified numerically with input parameters, such as mass, radius and magnetic field, 
characteristic to pulsars and magnetars. In Sec. 3, the frequency spectrum of axisymmetric 
global torsional nodeless vibrations under the action of the combined forces of shear elastic 
and magnetic field stresses is obtained and applied to ^-pole identification of the QPOs fre- 
quencies in flaring SGR 1806-20 and SGR 1900-1-14. The newly obtained results are briefly 
summarized in Sec.4. 
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2. Magnetic Alfven modes of global nodeless vibrations about etxis of 
magnetic dipole moment of a neutron star 

It is of common knowledge today that shortly before discovery of pulsars it has been re- 
alized that neutron stars must come into existence as compact sources of super strong dipolar 
magnetic field operating as a chief promoter of their electromagnetic activity. In particu- 
lar, it has been suggested that in view of the absence in neutron stars of thermonuclear 
source of electromagnetic emission from the star surface they could manifest their presence 
in the universe by the conversion of energy of hydromagnetic Alfven vibrations in a highly 
incompressible star mater into the energy of electromagnetic waves propagating out into 
the space (Hoyle, Narlikar & Wheeler 1964). In this work we investigate this proposition 
in the context of current development of the neutron star asteroseismology by proceeding 
from equations of magneto-solid-mechanics underlying study of non-compressional Alfven 
magnetic oscillations in quaking neutron stars. Such an approach implies that neutron star 
material, regarded as a perfectly conducting continuous medium threaded by a magnetic 
field frozen in the star on the final stage of gravitational collapse of its massive main se- 
quence (MS) progenitor, possesses mechanical properties of super dense strained solid highly 
robust to compressional distortions, rather than a flowing liquid, as is the case of MS stars 
whose theoretical asteroseismology relies of equations of fluid-mechanics and magnetohydro- 
dynamics. To this end, it worth emphasizing that contrary to the MS stars whose magnetic 
fields are generated by self-exciting dynamo process in a peripheral convective zone and 
the energy supply of this process comes from the central thermonuclear reactive zone (e.g. 
Parker 1979), in neutron stars there are no internal energy sources to drive relative motions 
of extremely dense conducting material which is under immense gravitational pressure so 
that mobility of neutron star matter is heavily suppressed. The common belief today is that 
magnetic fields of neutron stars are fossil (e.g. Fowlers & Rudcrman 1977, Bhattacharya & 
van den Heuvel 1991, Chanmugam 1994; Spruit 2008), that is, they have been inherited from 
a massive MS progenitors and amplified in the processes of implosive gravitational collapse 
proceeding under control of the magnetic fiux conservation. Also it is commonly agreed that 
external magnetic field of pulsars and magnetars is of dipolar symmetry, but the geometrical 
shape of the fossil magnetic magnetic field frozen-in the super dense motionless matter of 
the star remains fairly uncertain. The subject of present paper to formulate solid-mechanical 
variational scheme of computing frequency of Alfven oscillations in the neutron stars with 
arbitrary configuration of frozen-in fossil magnetostatic field. 

The equation of magneto-solid-mechanics describing coupled oscillations of material dis- 
placements u and magnetic field 5B about motionless stationary state of perfectly conducting 
elastically deformable neutron star matter in the presence of the frozen-in fossil magnetic 
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field is written as follows 

pii = ^[5jxB], 6} = ^[Vx6B], 5B = Vx[uxB], (1) 

pii = ^[V X [V X [u X B]] X B], V-u = 0, V • 5B = 0. (2) 

An 

Scalar multiplication of ([2]) by u and integration over the star volume leads to the equation 
of energy balance 

ly^rfV = i^|[V X [V X [u X B]] X B] . u dV. (3) 

This equation provides a basis for computing the frequency spectrum of nodeless vibrations 
by Rayleigh's energy variational method that has been used with success for analysis of 
shear nodeless vibrations driven by solely Hooke's elastic force (Bastrukov et al 2007; 2008). 
The key idea of this method is to use the following separable representation of the field of 
material displacements 

u(r,t) = a(r)a(t). (4) 

In a solid star undergoing global torsional nodeless oscillations about polar axis, the quadrupole 
and octupole overtones of which are pictured in Fig.l, is given by (Bastrukov et al, 2007a) 

u(r,t) = [u;{v,t) X r], uj = AtV/ Pe{cos9)a{t) (5) 

so that instantaneous displacements can be written as 

a, = A V X [rx(r)], = r'P,{cos9) : a,. = 0, a, = 0, = A,r\l - (2)1/2^^1 (6) 

where the arbitrary constant At which is eliminated from the boundary condition (Bastrukov 
et al, 2007a) 

u(r,t) = [u:{y,t) x r],=R = [II x R], O = AtW P,{cose)a{t) ^ At= [R'^'-%\ (7) 



^We remaind that equation for in line ([T]) follows from Maxwell equation for Faraday's induction 
5B = — cV X 5E where (5E = — (l/c)[(5v x B]. This link between fluctuating electric field (5E and constant 
magnetic field B frozen in the conducting flow follows from Ohm's law (Sj/cr = (5E + (l/c)[i5v x B] corrected 
for the extremely large, effectively infinite, electrical conductivity, cr — > cxd, a case of medium designated as 
a perfect conductor. Taking into account that (5v = li and eliminating time derivative from obtained in the 
above manner equation of the field-flow coupling (5B = V x [Jv x B] , one arrives at the last equation in the 
line (P). 
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On inserting (jlj) in ([3]) this latter equation is reduced to equation for a{t) having the form 
of equation of harmonic vibrations 

— = 0, ^ = ^— + 7Wa + /C^a = 0, (8) 

where the integral parameters of inertia M. and stiffness /C^ of magnetic Alfven vibrations 
are given by 

M=jpa^dV, ]Cm = ^ J[V x[\/ x[axB]]xB]-adV= ^ J[(B-V)a]'^dV. (9) 

In these last equations, the density p and frozen in the star magnetostatic fossil magnetic field 
B are considered to be intrinsic characteristics of equilibrium state of the star matter and 
input parameters of the method in question. The mass parameter Ai is positively defined 
quantity, whereas the sign of the stiffness /Cm depends on specific form of the magnetostatic 
fossil magnetic field B frozen-in the star. Thus, geometric configuration of the internal field 
is crucial to the question whether quake-induced shear perturbation resulting in fluctuations 
of differentially rotational displacements fllOl) are developed as oscillation mode, uj"^ > 0, or 
relaxation mode, uj"^ < 0, and within the above expounded method this issue, that has been 
a subject of controversy (Levin 2006, Glampedakis et al 2006; Watts & Strohmayer 2007), 
cannot be resolved with no computing Km for each imaginable form of B. 

The large scale external magnetic fields of radiative magnetospheres of pulsars and 
magnetars are commonly thought of as produced by the magnetic dipole moment of under- 
lying neutron star and it is generally believed that internal magnetic field of the neutron 
star has a strong poloidal dipolar component. With this in mind and for the reason of 
computational feasibility, in the remainder of this section we examine the above variational 
approach by considering an admittedly idealized neutron star model with uniform poloidal 
internal magnetostatic magnetic field B directed along the polar axis z undergoing global 
torsional oscillations about dipolar magnetic moment axis which are insensitive to the above 
compositional, core-crust, stratification of the neutron star. The explicit form of spherical 
components of such a field inside the star are given by 

Br = BC, Be = -B{l-Cy^\ 5^ = 0, r < R, ( = cos9 (10) 

and components of dipolar configuration outside the star reads 

5. = 5(^)'C, Be = {1 - ey^', B, = r>R (11) 

where R is the the star radius (e.g. Chandarsekhar & Fermi 1953). Computation of integrals 
for Ai and /Cm yields (Bastrukov et al, 1997) 
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and for the angular frequency of toroidal magnetic Alfven mode of global torsional nodeless 
oscillations we obtain 



a;(o4) = loa 



2 . . 2^ + 3 



- 1) 



1/2 



21 -l_ 

The fundamental frequency of Alfven oscillations ua can be conveniently written in the form 



where M and i? stand for the mass and the radius of star. It may be worth noting that 
normal component of magnetic field under consideration is continuous on the surface while 
the tangential component remains discontinuous and, thus, admits, in accord with standard 
boundary condition of electrodynamics, the surface current, provided that on the surface 
there is an excess of likely-charged particles distributed with the surface charge density a. 
The considered magnetic field inside the star is identical to that produced by uniformly 
charged spherical shell of radius R set in the rotation about polar axis with constant angular 
velocity Vl and the absolute value of this field is given hy B = {2/3)ailR (Griffiths 1981). 
So, the above inferences regarding frequency of nodeless Alfven oscillations in the star vol- 
ume are not affected when electrodynamic conditions inside the neutron star model under 
consideration are compatible with those for the latter case of uniformly charged spherical 
shell. 

In Fig. 2, frequencies (left) and periods (right) of the toroidal Alfven mode are plotted 
as functions of multipolc degree £ for a solid star with parameters typical for pulsars and 
magnetars clearly showing general trends of u (in Hz) and P (in seconds) as a function of 
multipole degree of oscillations and order of their magnitude. 

In the above computations we have used toroidal (axial) vector field of differentially- 
rotational nodeless material displacements which is one of two fundamental solutions of the 
vector Laplace equation 

V^a(r) = 0, V • a(r) = (15) 

built on fundamental solution of the scalar Laplace equation 

V\(r) = 0, x{r)^r'Pe{C), C = cos^. (16) 

The second fundamental solution is given by the even-parity poloidal (polar) vector field 
(Bastrukov et al 2007a) 

ap = X V X [rx(r)] = A^Vx(r), xW = r^P^(C) (17) 
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which is irrotational: V x = 0. The integral parameters of inertia Ai and stiffness JCm in 
the poloidal mode are given by 

Mioa^,) = ^^pKR"'"-'^^ ^-(o«^) = B'AlR^'-'^^^^ (18) 
and for the frequency spectrum of poloidal magnetic Alfven mode we get 

nx o ^ / ^ .2£-\-l r, V] B^R 

.(,4) = .w-i)—, .j = ^ = _ = _ (19) 

The spectral formulae like above obtained for toroidal and poloidal Alfven mode^ are central 
to theoretical asteroseimology of pulsars and magnetars in the sense they provide the basis 
for interpretation of observable quasi-periodic oscillations as produced by quake-induced 
vibrations restored by Lorentz force. 

For the purpose of our further consideration we remind that the Lorentz force density, 
f = (l/c)[j X B], with the Ampere's current density, j = (c/47r)V x B, having the form 



— [[V X Bl X Bl = — 

ATX^^ ^ ^ 4:71 



(B • V)B - ^VB^ 



(20) 



can be represented in terms of Maxwell tensor of magnetic field stresses as (e.g., Mestel 1999) 

fi = ^kTik, = — [Bi Bk + Bk Bi — {Bj Bj)6ik] ■ (21) 
ovr 

In an extremely dense matter of solid stars of the finale stage the effects of magnetic buoyancy 
are heavily suppressed and, hence, in the motionless matter of the stationary state of the solid 
star with the above frozen-in uniform magnetic field fi = ff{B = const) = 0. However, the 
Lorentz force surely comes into play as a result of quake induced perturbation which leads to 
coupled fluctuations of material displacements u and magnetic field (5B in accordance with 
equation for field- flow coupling 5B(r, t) = V x [u x B]. With this in mind, the equation of 
dynamics ([2]) can be represented in the equivalent tensor form, to wit, in terms of tensor of 
fluctuating magnetic field stresses (e.g.. Franco, Link & Epstein 2000) 

pui = Vk Tik, r,k = -^[B, 5Bk + Bk 5B, - {B, 5B,)5,k], 5B = V x [u x B] (22) 



■^To the best of our knowledge, classification of Alfven vibrational modes as poloidal and toroidal has 
been introduced by Chandrasekhar (1956) in the context of hydromagnetic oscillations of fluid sphere. Gen- 
eral properties of poloidal and toroidal fields are extensively discussed in (Chandrasekhar 1961; Ferraro & 
Plumpton 1961). 
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where B is the in-advance-given stationary magnetic field. In the text section we show that 
above obtained spectral equation f|T3|) for the toroidal Alfven mode of global differentially 
rotational oscillations of the star matter about homogeneous field frozen-in the star can be 
regained on the basis of this last representation of the Lorentz force density. 

In Fig. 3 we plot the frequency and period of Alfvenic modes uj{oai) upon multipole 
degree i of nodeless oscillations of material displacements in a neutron star whose internal 
magnetic field might be approximated by axisymmetric poloidal field of uniform shape. The 
presented absolute values for the frequencies show that they fall in the realm of observable 
QPOs in fiaring SGR's 1806-20 and 1900-1-14 and, thus, suggesting that Lorentz restoring 
force must be taken into account when studying global seismic vibrations of magnetars. 
Together with this it worth emphasizing that obtained one-parametric spectral formula for 
toroidal Alfven vibration mode in model under consideration does not match the general 
trends in observed QPOs frequency as a function of multipole degree of torsional vibrations 
(Watts & Strohmayer 2007) because the slope of computed frequency uj{oa\) as a function 
of i is different from the slope of the overall trends of the observed QPOs frequencies. In the 
next section, we approach to this problem by considering in some details the global nodeless 
torsional oscillations of a neutron star model with the above uniform internal magnetic field 
under the action of combined forces of shear elastic and magnetic field stresses. 

3. Torsional nodeless vibrations under the joint action of bulk forces of elastic 

and magnetic field stresses 

The equation of Newtonian magneto-solid-dynamics appropriate for non-compressional 
{6p = — pV ■ u = 0) shear vibrations of an elastic medium of infinite electrical conductivity 
are 

pili = Vfe (Tik + Vfc Tik. (23) 

The first term in the right part of fl23l) is the bulk force of elastic stresses obeying the Hooke's 
law 

aik = 2/iUifc, Uik = ^(ViMfc + VkUi), Ukk = VkUk = (24) 

where /i is the shear modulus of stellar matter linearly relating quake-induced shear stresses 
and resulting shear deformations or strains Uik- The second term is the above defined 
Lorentz force represented in terms of fiuctuating magnetic field stress Tik- Similar equations 
have recently been considered in works (Piro 2005; Glampedakis, Samuelsson & Andersson 
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2006; Lee 2008). The conservation of energy is controlled by equation 
d f pix^ixi f 1 

^ / -^-^ dV = - [cTifc + Tik] iiik dV, iiik = -[Viiik + VkUi] (25) 

which is obtained after scalar multiplication of ( l23l) by iii and integration over the star vol- 
ume. To compute eigenfrequency of torsional nodeless oscillations we again take advantage of 
the Rayleigh's energy method. On inserting separable representation of the field of material 
displacement 

Ui{r,t) = ai{r)a{t) (26) 

in equation of energy balance we again obtain equation of harmonic oscillations of temporal 
amplitude a{t): 

^ = 0, n = ^ + ^, ^ >[«(t) + /C«(t) = 0, /C = /Ce + /C„ (27) 

Analytic form for the inertial Ai and for the stiffness JCe of shear elastic oscillations derived 
in (Bastrukov et al 2007a; 2007b) are given by 

M = j pai{r) ai(r)rfV, Ke = ^ j n[Vkai + V^ak] [Vfcfli + V^afc] dV (28) 

and for the stiffness of magnetic Alfven shear oscillations we get 

= ^/{a|V,Wa.-B.a,)]+B.|V,(B,a.-B.a,)l}(V.a, + V.a.)<iV 

It should be noted that this last equation can be derived from the standard MHD equations 
(Chandasekhar 1961). Computation of integrals with Bi given by ([3]) and toroidal field 
defined by ([9]) yields 

M,M ^ 4.M?«-» p^flyl^. ^ (30) 

K,(„e;) = 4xM?fi^'«|^. A„(oa')^BM?fl^'^' g+j><^^^:^^'^ (31) 

and for the measurable in Hertz total frequency u = uj/2'k (with uj = y^lC/Ai) we obtain 

z/(ot£) = [^e(oeD + z/^(o4)]'/', (32) 
z/e(oe*) = z/e[(2£ + 3)(£-l)]l/^ z/, = ^, u;, = ^, Ct = J^, (33) 



2 , ^ 2£ + 3 



ii' - 1) 



2^-1 



1/2 



(^A VA B 
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In Fig.4, the fractional frequency of elastic i>e{oG*t)/^e oscillations as a function of multipole 
degree £ is plotted in juxtaposition with fractional frequency VmioCi^i) / i^A of magnetic Alfven 
oscillations. One sees that the lowest overtones of both elastic and magnetic Alfven modes are 
of quadrupole degree, 1 = 2. At £ = 1, both parameters of elastic and magneto- mechanical 
rigidity cancel, /Ce(oei) = and /Cm,(ofli) — ^"^^ ^^e mass parameter equal to the moment 
of inertia of rigid sphere, }A = {2/b)MR^ . It follows from Hamiltonian that in this dipole 
case a star sets in rigid body rotation, rather than vibrations, about axis of magnetic dipole 
moment. 



The obtained spectral equation for total frequency can be conveniently represented in 
the following form 



v{oti) = i/e[(2£ + 3)(£-l)]i/2 



l + (3 



21-1 



1/2 



(35) 



where 



c4 ^ ^ 4 

P .9 -.7 ,2 



(36) 



In Fig. 5 we plot the last two-parametric spectral equation for the total frequency vloti) 
as a function of multipole degree i of torsion nodeless vibrations computed with indicated 
values of parameters Ve a-nd f3 carrying information about mechanical and electrodynamical 
properties of the neutron star matter which are adjusted so as to reproduce the observable 
frequency of QPOs (symbols) during the flare of SGRs 1806-20 and 1900+14; the data from 
(Watts & Strohmayer 2007, Samuelsson & Andersson 2007). This figure demonstrates that 
the detected QPOs can be consistently explained from the viewpoint of the the considered 
model as produced by global torsional nodeless vibrations when and only when combined 
forces of shear elastic and magnetic field stresses come into play in a coherent fashion. On 
the other hand, in our previous study reported in recent paper (Bastrukov et al, 2008a) 
it has been shown that this set of QPOs data can be properly described, with the same 
degree of accuracy, on the basis of two-parametric spectral formula that has been derived 
on the basis of a two-component, core-crust, model of quaking neutron star presuming that 
detected QPOs arc produced by axisymmetric torsional nodeless seismic vibrations driven 
by a solely clastic restoring force and locked in the peripheral finite-depth scismogenic layer. 
The truth is, most probably, somewhere in between and in order to attain more definite 
conclusions further investigations, both theoretical and observational, are needed. Prom a 
computational argument, the nodeless torsional oscillations entrapped in the neutron star 
crust as well as in the star models with non-uniform axisymmetric internal magnetic field 
and non-homogeneous profile of shear modulus, requires a more elaborate mathematical 
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treatment. To avoid destructing attention from the newly obtained results presented here 
[some of which are of interest, as is hoped, for general theoretical seismology (Lay & Wallace 
1995; Aki & Richards 2002)], we postpone a discussion of these latter cases to a forthcoming 
article. 

A special comment should be made regarding the link between neutron stars and atomic 
nuclei which can be considered to be similar objects as far as mechanical properties of degen- 
erate nucleon material of normal nuclear density are concerned. It follows from the nuclear 
solid-globe model of giant resonances (within the framework of which these fundamental 
modes of nuclear excitations arc properly described in term of shear clastic vibrations of an 
ultra small spherical piece of nucleon Fermi-solid (regarded as continuous matter) that the 
shear modulus is given by ~ 10^^ dyn cm~^ (e.g., Bastrukov et al, 2008b). Making use of 
this value of /i in the expression for parameter /3 = [B'^/Anij] ~ 0.5 — 0.6 entering the above 
derived spectral equation for QPOs frequency, we can get an independent estimate for the 
intensity of the magnetic field frozen in the star. With the above values of f3 and /j, one finds 
that B falls in the range 10^^ < B < 10^^ Gauss, that is, in the realm typical of magnetar 
magnetic fields. This latter inference demonstrates the theoretical potential of the astero- 
seimology showing how physical interpretation of the QPOs as produced by quake induced 
torsional vibrations of neutron star can be used to extract information about properties of 
the neutron star matter. 



4. Summciry 

There is a common belief today that gross features of the asteroseimology of pulsars 
and magnetars can be understood on the basis of a solid star model presuming that quake- 
induced shear vibrations restored by bulk forces of intrinsic stresses of different physical 
nature are governed by described by solid mechanics or elastodynamics (e.g., Hansen & van 
Horn 1979; McDermott, van Horn & Hansen 1988; Bastrukov, Weber & Podgainy 1999; 
Bastrukov et al 2007a). This point of view is quite different from theoretical approach to the 
asteroseismology of the main sequence stars at the base of which lies the liquid star model 
whose vibrations are treated within the framework of fiuid-mechanical theory of continuous 
media, as is the case of helioseismology (e.g., Christensen-Dalsgaard 2002). 

The main purpose of this work was to examine the magneto-solid-mechanical varia- 
tional method of the asteroseismology of neutron star by probing its interior with non-radial 
global differentially rotational, torsional, vibrations with nodclcss toroidal field of material 
displacements, which are insensitive to compositional stratification of the star matter. Bear- 
ing in mind that external magnetic fields of pulsars and magnetars are commonly thought 
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of as produced by magnetic dipole moment of underlying neutron star we have considered a 
model of a neutron star with perhaps simplest, from the viewpoint of computational feasi- 
bility, imaginable configuration of the magnctostatic fossil magnetic field, pictured in Fig.l. 
Proceeding from this admittedly idealized model, the two-parametric spectral equation for 
the frequency of global torsional vibrations has been derived in analytic form showing the 
larger multipole degree I of torsional nodeless vibrations the higher is the frequency uii^oti). 
The application of the obtained two-parametric spectral formula to modal analysis of QPOs 
during the flare of SGR 1806-20 and SGR 1900+14 shows that data on the QPOs frequen- 
cies with I from interval 2 < ^ < 20 can be consistently interpreted as produced by global 
torsional nodeless vibrations restored by combined forces of shear elastic and magnetic field 
stresses. This inference is, of course, suggestive rather than conclusive, in view of adopted a 
highly idealized configuration of fossil magnctostatic field frozen in the star, and too much 
remains to be done to be at all certain of interpretations suggested for QPOs in the X-ray 
flux during the giant flares of the above magnetars. 
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Fig. 1. — Fiducial model of a solid star with uniform magnetic field inside and dipole outside 
undergoing global torsional nodeless oscillations in quadruple and octupole overtones. 
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Fig. 2. — Frequencies (left) and periods (right) of toroidal magnetic Alfven vibrational mode as 
functions of multipole degree i in the neutron star models with pointed out intensities of uniform 
magnetic field inside and dipole outside. 
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Fig. 3. — Computed frequencies and periods of poloidal and toroidal Alfven modes in the neu- 
tron star models undergoing global spheroidal and torsional shear nodeless vibrations, respectively, 
driven by Lorentz restoring force. 




Fig. 4. — Fractional frequencies v/vq of toroidal Alfven mode (with vq = va) and torsional elastic 
mode (with fo = t'e) in a solid star undergoing differentially rotational nodeless oscillations about 
polar axis of uniform inside and dipole outside magnetic field. 
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Fig. 5. — Theoretical fit of data (symbols) on QPOs frequencies during the flare of SGRs 1806-20 
and 1900+14 by the above spectral equation for the frequency of axisymmetric torsional nodeless 
vibrations about axis of dipole magnetic moment of magnetar under the joint action of bulk forces 
of elastic and magnetic stresses. 



